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ANALYTICAL CHEMICAL KINETIC STUDY 
OF THE EFFECT OF CARBON DIOXIDE AND WATER VAPOR 
ON HYDROGEN-AIR CONSTANT-PRESSURE COMBUSTION 

By Wayne D. Erickson and George F. Klich 
Langley Research Center 

SUMMARY 

Numerical solutions have been obtained for the finite-rate constant -pressure com- 
bustion of stoichiometric hydrogen-air mixtures in the presence of small to moderate 
amounts of carbon dioxide (CO2) and water vapor (H2O). Computations have been carried 
out for initial mixture temperatures of 1150 K, 1250 K, and 1500 K at a pressure of 
1 atm, with additional computations for pressures of 0.5 atm and 2 atm at an initial mix- 
ture temperature of 1250 K. The results show that the influence of CO2 and H2O on the 
induction time, for the condition of no free radicals (O, H, and OH) initially present, is 
quite small (an increase up to about 16 percent) for initial temperatures of 1250 K or 
higher at pressures of 1 atm or less. For an initial mixture temperature of 1150 K at 
1 atm or 1250 K at 2 atm, the induction time without initial free radicals is increased 
considerably (nearly 60 percent) by the presence of moderate quantities of CO2 and H2O. 
These increases in induction times due to the presence of initial amounts of CO2 and H2O 
are offset somewhat when initial small amounts of free radicals (O, H, and OH) are 
included in the calculations. In some cases the inclusion of initial quantities of free 
radicals causes a reduction in induction time due to initial amounts of CO2 and H2O. The 
reaction time was measurably increased by the presence of CO2 and H2O for all test 
cases, but a large fraction of this increased time can be accounted for by the reduction 
in temperature due to the added CO2 and H2O. This study suggests that although all the 
conditions for hydrogen-burning hypersonic ramjet engine tests in a combustion-heated 
wind tunnel cannot be matched to clean-air flight conditions, the chemical kinetic effects 
of CO2 and H2O are small enough to allow useful interpretation of test results for initial 
temperatures of 1250 K or higher and pressures near 1 atm or less. 

INTRODUCTION 

A proposal has been made to carry out hydrogen-burning hypersonic ramjet engine 
tests in a combustion-heated wind tunnel in which the test stream is composed of the com- 
bustion products formed by burning a hydrocarbon fuel in oxygen- enriched air. The 



Langley 8-foot high-temperature structures tunnel, if modified to include an oxygen- 
enrichment system, would be a facility of this kind. The use of such a simulated test 
stream raises the question of the influence of the tunnel-heater combustion products, 
mainly carbon dioxide (CO2) and water vapor (H2O), on the chemical kinetic processes 
in the test engine. 

There have been, of course, many analytical and numerical studies of hydrogen- 
oxygen (H2-O2) combustion processes. Some of the more recent are listed as refer- 
ences 1 to 5. The effects of vitiation have also been the subject of studies such as ref- 
erences 6 to 9. It may be noted, however, that these studies have emphasized the ignition 
phase of the combustion process and do not generally consider the overall process which 
includes that final approach toward equilibrium. The gas temperature range considered 
in these previous studies are also generally higher than would apply to a tunnel such as 
the Langley 8-foot high-temperature structures tunnel, which is heated by methane 
combustion. 

The purpose of this paper is to examine the influence of small to moderate amounts 
of CO2 and H2O on the induction time and reaction time for hydrogen-air combustion 
under constant pressure. The amounts of CO2 and H2O chosen for these calculations are 
representative of the concentrations of these two species in a tunnel test stream gen- 
erated by oxygen-enriched methane-air combustion. For comparison, computations are 
also made for cases with no CO2 and H2O initially present. The calculations are based 
on the assumption that the hydrogen and oxidizer streams are completely mixed at the 
entrance of the combustor. 


SYMBOLS 

Aj preexponential reaction rate coefficient 

ai j relative third-body efficiency of species i in reaction j 

Bj power on T in equation (8) 

ci mass fraction of species i 

Ej activation energy of reaction j 

h static enthalpy (including enthalpy of formation) 

hi static enthalpy of species i 
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static enthalpy at entrance of combustor 


P 

R 

T 

T e 

T m 

T 0 

t 


u 

Wi 

Xi 

x 



number of chemical species and finite-rate reactions, respectively 

specific chemical species and finite-rate reaction, respectively 

equilibrium constant in terms of molar concentrations for reaction j 

reaction rate constant for reaction j 

reaction rate constant for reverse of reaction j 

third-body species 

specific third-body species 

pressure (1 atm = 101.325 kN/m2) 

universal gas constant 

temperature 

equilibrium adiabatic flame temperature 
mean temperature, — (t o + T e) 
initial temperature at entrance of combustor 
time (t = 0 at entrance of combustor) 
velocity 

molecular weight of species i 
general symbol for chemical species i 
distance along combustor, measured from entrance 
fractional temperature rise (see eq. (12)) 
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T 


At 


T r 


i,] 



v i+i,i 


P 


induction time without initial quantities of free radicals O, H, and OH 
(time at which 6 = 0.05) 

difference in induction time with and without free radicals 

reaction time (the total reaction time to 6 = 0.95 less the induction time 
at which 6 = 0.05) 

stoichiometric coefficient of reactant species i in reaction j 

stoichiometric coefficient of product species i in reaction j 

index in equation (9) which is unity for all three -body reactions and 
zero for all bimolecular reactions 

gas density 


ANALYSIS 


The assumptions and equations for computing the supersonic flow properties along 
a constant-pressure combustion chamber are outlined in this section. It is assumed that 
the flow is quasi-one -dimensional, inviscid, and adiabatic. The hydrogen and test stream 
are assumed to be completely mixed just prior to the beginning of the combustion process. 


Thermodynamic Relations 

On the basis of the foregoing assumptions, the conservation of momentum 


and the conservation of energy 


du , dp n 
pu — + = 0 

dx dx 


( 1 ) 



( 2 ) 


combine to show that , when ^ = 0 

’ dx 

bustor is constant — that is, 


in equation (1), the static enthalpy along the com- 

^ = 0 (3) 

dx 


It also follows that the flow velocity along the combustor is constant - that is, 

*i=0 

dx 


( 4 ) 
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and therefore the dwell time t of a given mass of reacting gas mixture is directly pro- 
portional to x, the distance from the entrance of the combustor. Since the static 
enthalpy of the reacting gas mixture, composed of I chemical species, in the combustor 
is independent of time, one equation which must be satisfied at every point in the com- 
bustor is 

I 

£ Cihi(T) = h D (5) 

i=l 


where h Q is the static enthalpy at the entrance of the combustor. Equation (5) applies 
to any point in the combustor and indicates how the temperature relates to a particular 
chemical composition expressed by a set of mass fractions. The expression which 
relates the gas density to the temperature and chemical composition at a particular pres- 
sure is the equation of state 



( 6 ) 


which must be solved along with equation (5). 

For finite-rate combustion the differential equations which represent the various 
chemical reactions must be solved simultaneously in a stepwise manner along with equa- 
tions (5) and (6). These differential rate equations are discussed in the next section. 


Chemical Kinetic Relations 


The reacting gas system is assumed to contain J finite-rate reactions, all of 
which may be expressed as 

1 k 1 

I t . Z " i ,)* 1 (i = 1 ’ 2 J) (7) 

“3 i=l 


i=l 


where the rate constant kj for the jth reaction is represented by 

kj = AjT B i 


(i = 1, 2, . . . J) 


( 8 ) 


The rate of change of the mass fraction of the ith chemical species in the jth reac- 
tion, represented by equation (7), can be expressed as 

1 I- a, .cV^M 




(J - 1, 2, ... J) 


(9) 


The last term in this equation accounts for the homogeneous catalysts, say M, in three- 
body reactions. The index a^j is the relative third-body efficiency of species i in 
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reaction j if the ith species is to be included as a third-body catalyst in the jth reac- 


tion, but is zero if that species is to be excluded. The quantity u 


1+1,3 


is also an index 


and is equal to unity for all three-body reactions. For bimolecular reactions, j 

is zero with all values of a^j set equal to unity. 

The total rate of change of the mass fraction of the ith species in the reacting gas 
mixture is obtained by summation over all J reactions; that is, 


dc 

dt 


1 = V/ dc i 


3=1 J 


(i = 1, 2, ... I) (10) 


where the summation includes both forward and reverse reactions. The set of equations 
to be solved includes the I differential equations indicated by equation (10), along with 
equations (5), (6), (8), and (9). The unknowns are reduced to Cj, Cg, . . . Cj and T, 
after equations (6) and (8) are used to eliminate p and kj. The initial values of the 
mass fractions and temperature are specified at the entrance of the combustor so that a 
forward integration of this set of equations can be carried out along the length of the com- 
bustor by use of a modified Runge-Kutta scheme. 


Chemical Kinetic Scheme 

The chemical kinetic scheme used in this investigation included the following 
reactions: 


K 1 

(1) OH + H 2 Z H 2 0 + H 
k -l 
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(6) O + O + M — 0 2 + M 


(2) H + 0 2 Z 2 OH + O 
k -2 


k 7 

(7) 0+H + M=' OH + M 
k -7 


k 3 

(3) O + H 2 r OH + H 
k -3 


(8) H + OH + M r 8 H 2 0 + M 
k -8 


k 4 

(4) OH + OH r h 2 o + o 
k_ 4 


k 9 

(9) CO + OH - C0 2 + H 

k -9 


(5) H + H + M r 5 H 2 + M 
k -5 


k in 

(10) H + 0 2 + M Z H0 2 + M 
k -10 


6 



(11) H0 2 +H 2 = n H 2 0 2 +H (13) H 2 +0 2 r 13 OH + OH 

k_n k_ 13 

k-i o 

(12) H 2 0 2 + M - OH + OH + M 

k -12 

The values of the reaction rate constant kj (for j = 1, 2, . . . 13) used in this study 
are listed in table 1 along with the sources (refs. 10 to 18). The reverse reaction rate 
constants, not shown in table 1, were calculated from the equation 

k -j = k*: (11) 

where the equilibrium constant for the jth reaction Kj was derived from thermochemi- 
cal data of reference 19. For the three-body reactions the values of kj are indicated 
in table 1 for a specific species, say M’ , and the relative third-body efficiencies of the 
remaining species are shown in table 2, along with the sources (refs. 15, 17, and 20 to 24). 

CALCULATIONS 

Three test-gas mixture compositions were chosen in order to determine the influ- 
ence of moderate amounts of C0 2 and H 2 0 on the rate of H 2 -air combustion. Table 3 
lists these three combustion-chamber-inlet mixtures which include added H 2 . All mix- 
tures are stoichiometric with respect to H 2 and the amount of 0 2 to form H 2 0. The 
nitrogen (N 2 ) is assumed chemically inert. Mixture R is a reference mixture which 
contains no C0 2 or H 2 0. The relative amounts of C0 2 and H 2 0 in mixtures A and B are 
based on complete stoichiometric combustion of CH4 and air (in the heater section of the 
test facility) with 0 2 added to bring the stream 0 2 mole fraction to 0.2 before H 2 addi- 
tion. The species included in mixtures A and B simulate products of combustion with 
tunnel stagnation temperatures of 1390 K and 2200 K, respectively, which are approxi- 
mately the lower and upper limits for methane-air combustion-heated tunnels. 

The three mixtures were used as the basis for calculating the finite-rate process 
along a combustion chamber for various combinations of initial temperature and pres- 
sure. The cases calculated are listed in table 4. For each of these 15 test cases the 
equilibrium adiabatic flame temperature was also computed. The finite-rate calcula- 
tions were allowed to proceed in time for each case until at least 95 percent of the maxi- 
mum possible temperature rise (the equilibrium adiabatic flame temperature less the 
initial temperature) was achieved. The foregoing set of calculations were first carried 
out with the assumption of no initial free radicals (O, H, and OH). The same set of 
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calculations were then performed with the assumption of an equilibrium amount of free 
radicals in the separate hydrogen- and oxygen- containing streams at the initial tempera- 
ture To and pressure p in the combustor prior to mixing. The sums of the free 
radicals in the two streams were added together to give the initial mole fractions of free 
radicals at the entrance of the combustor. These quantities are listed in table 5. The 
results of these calculations without and with initial free radicals are shown and dis- 
cussed in the next section. 


RESULTS AND DISCUSSION 


Calculations have been made for the 15 cases indicated in table 4. Figures 1 to 15 
show the time variation of the species mass fraction. In addition, the equilibrium condi- 
tion in terms of mass fractions is also shown on the far right of these figures. 


Figure 1 shows the detailed behavior of each reactant and product species for 
case 1, which was computed for an initial temperature of 1250 K, a pressure of 1 atm, 
and mixture R containing no initial CO 2 or H 2 O and no initial free radicals. The attain- 
ment of maximum levels and subsequent decay toward equilibrium is observed for the 
intermediate species (OH, O, H, HO 2 , and H 2 O 2 ), while H 2 O grows at the expense of 
depleting O 2 and H 2 . The level of N 2 is fixed since it is inert. Figures 2 and 3 are based 
on the same initial temperature and pressure, but are for successively larger initial 
amounts of CO 2 and H 2 O. These two figures include the variation of CO 2 and CO and 
indicate an initial concentration of H 2 O. 

Figures 4 to 6 represent similar results but for a pressure of 0.5 atm, with corre- 
spondingly longer times, whereas figures 7 to 9 are for a pressure of 2 atm. The reac- 
tion times for the 2-atm cases are, as expected, considerably shorter than for the lower 
pressure cases. Figures 10 to 15 show the effect of different initial temperatures 
(1500 K and 1150 K) on the species time histories. 


In order to obtain a quantitative assessment of the influence of the addition of CO 2 
and H 2 O on the induction time and reaction time, the degree of reaction is defined by the 
equation 


6 = 


T - T 0 
T e -T 0 


( 12 ) 


where 6 ranges from zero at the entrance of the combustor toward unity as equilib- 
rium is approached. The time variations of the fractional temperature rise 6 for the 
15 cases are shown in figures 16 and 17. 

In this study the induction time is defined as the time at which 6 = 0.05, and the 
total reaction time is defined by 6 = 0.95. The induction times calculated for the 
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reference mixture R for pressures of 0.5 and 1 atm have been compared with experi- 
mental and calculated results presented in references 1 and 25. It was found that the 
present calculations agree within the range of uncertainty of the reported induction times. 

Table 6 lists the calculated induction times for various initial temperatures, pres- 
sures, and amounts of CO2 and H2O when (1) there are no initial free radicals (O, H, and 
OH) and when (2) the initial free radical concentration is based on the sum of the equilib- 
rium amounts in the two separate streams (hydrogen- and oxygen-containing streams) 
prior to mixing. It is noted from table 6 that, with no initial free radicals, the addition 
of CO2 and H2O causes an increase in induction time of no more than 16 percent for ini- 
tial temperatures of 1250 K or higher at pressures of 1 atm or less. However, as the 
pressure is increased to 2 atm at 1250 K, it is seen that the induction time increases by 
almost 60 percent for mixture B, and an even larger increase in induction time due to 
CO2 and H2O addition is noted for the 1150-K, 1-atm case (without initial free radicals). 

An attempt to indicate the influence of an initial concentration of free radicals was 
made by calculating the induction times for the initial free radical concentrations listed 
in table 5. These results are presented in table 6 along with the results for no initial 
free radicals. It is noted that the increases in induction time in all cases just discussed 
are offset to some degree by an initial free radical concentration. In some cases the 
effect of an initial free radical concentration causes the induction time to be less for the 
mixtures containing initial quantities of CO2 and H2O than for the reference mixtures. 

In an actual test, it is difficult to determine the amount of free radicals initially 
present. Figure 18 was therefore prepared in order to correlate the relative reduction 
in induction time - — as a function of the total mole fraction of initial free radicals. 

T 

The sums of the mole fractions of O, H, and OH are used because the growth for each free- 
radical species is strongly coupled to the other free radicals. It can be seen from fig- 
ure 18 that the correlation of as a function of the sum of the initial mole frac- 

T 

tions O, H, and OH is very good for T 0 = 1250 K where a number of results have been 
obtained. The results for T 0 * 1150 K and T 0 = 1500 K are rather limited but sug- 
gest a reasonable correlation also. 

The total reaction time required to achieve a 95-percent temperature rise (0 = 0.95) 
less the induction time (at which 6 = 0.05) for each of the 15 cases is given in table 7(a) 
as reaction time T r . Calculations with and without initial free radicals show that this 
reaction time is independent of the initial free radical concentration. The reaction time 
is, however, influenced by the presence of initial quantities of CO2 and H2O. An examina- 
tion of the calculated equilibrium temperatures as shown in table 7(b) suggests that the 
increases in T r due to initial amounts of CO2 and H2O may result to a large measure 
from the reduction in an average reaction temperature. Figure 19 shows the reaction 
time including a pressure-dependent term T r p*'® as a function of the reciprocal of a 
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mean temperature l/T m . It can be seen that the expression log r r p^*® = + 0.040 

correlates the results to within ±10 percent. This figure does suggest that a large part 
of the increase in r r is due to the reduction in the mean reaction temperature T m . 

CONCLUDING REMARKS 

Numerical solutions have been obtained which show the chemical kinetic effect of 
small to moderate amounts of carbon dioxide (CO 2 ) and water vapor (H 2 O) on hydrogen- 
air combustion under constant pressure. These results were obtained for mixtures 
which cover a range of temperature, pressure, and chemical composition in an oxygen- 
enriched methane-heated wind tunnel which might be used for a hydrogen-burning ramjet 
engine test. 

The results of this study for initially uniformly mixed systems show that the induc- 
tion time is increased by no more than 16 percent by the presence of moderate amounts 
of CO 2 and H 2 O for initial temperatures of 1250 K or higher at pressures of 1 atm or less 
when no free radicals are initially present. The initial presence of CO 2 and H 2 O at 
lower initial temperatures or higher pressures causes rather large increases in induc- 
tion time for no initial free radicals. However, the initial presence of rather small 
quantities of free radicals offsets to some degree these increases in induction times. In 
some cases an allowance for the initial presence of free radicals causes a net reduction 
in induction time for mixtures containing CO 2 and H 2 O relative to the reference mixtures. 
The increase in reaction time (total reaction time less induction time) due to the pres- 
ence of CO 2 and H 2 O appears to result primarily from the diluting effect which reduces 
the mean reaction temperature. 

On the basis of this limited set of calculations, it appears, from a purely chemical 
kinetic point of view, that results from hydrogen-burning ramjet engine tests in an 
oxygen-enriched methane -heated wind tunnel would not be grossly different from clean- 
air tests provided the combustion inlet temperature is 1250 K or higher and the pressure 
is 1 atm or less. Additional calculations, however, would be required to show more pre- 
cisely the temperature and pressure limitations and allow interpretation of engine test 
results. The influence of the mixing process, which was not studied herein, must of 
course also be examined. Another point to be noted is that the values of the reaction rate 
constants used in the present kinetic scheme are of course subject to improved accuracy, 
but it is expected that the trends and approximate level of the results that are presented 
are correct. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 10, 1970. 
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TABLE 1.- REACTION RATE CONSTANTS 


i 

Reaction 

Reaction rate constant 
(a) 

Source 

l 

k i 

OH + H 2 = H 2 0 + H 

kj = 2.3 X 10 13 exp(-5200/RT) 

Dixon-Lewis, Wilson, and 
Westenberg (ref. 10) 

2 

k 2 

H + Oo OH + O 

2 k-2 

k 2 = 1.0 X 10 14 exp(-16 000/RT) 

Recommendation by 
Richard S. Brokaw, 
NASA Lewis Research 
Center (1967) 

3 

k 3 

O + H 2 r OH + H 
2 k -3 

k 3 = 4.3 X 10 13 exp(-10 200/RT) 

Wong and Potter (ref. 11) 

4 

k 4 

OH + OH = HoO + O 
k -4 2 

k 4 = 7.6 X 10 12 exp(-1000/RT) 

Kaufman and Del Greco 
(ref. 12) 

b 5 

kg 

H + H + M H 9 + M 

k-5 2 

k -5,Ar = 1 - 18 x 1q12 t °’ 5 ° exp(-94 500/RT) 

Myerson, Watt, and Joseph 
(ref. 13) 

b 6 

k 6 

o + o + Mr o 2 + m 
k -6 

k R 0 = 8.2 X 10 18 T" 1,22 
6 ,°2 

Kiefer and Lutz (ref. 14) 

b 7 

k 7 

O+H + Mr OH + M 
k -7 

k 7j M = 6.°Xl° 14 

Kaskan and Browne (ref. 15) 

b 8 

k 8 

h + oh + m: h 2 o + m 

k -8 

k 8,H 2 0 = >•* X 1022 T ' 1,5 

Kaskan and Browne (ref. 15) 

9 

k 9 

CO + OH ^ C0 2 + H 

-9 

kg = 3.1 X 10 11 exp(-600/RT) 

Dixon-Lewis, Wilson, and 
Westenberg (ref. 10) 

b 10 

H + o 2 + m r 10 ho 2 + M 
k -10 

k 10 Ar = 8,6 x 1q14 exp(1280/RT) 

Getzinger and Schott (ref. 16) 

11 

k ll 

ho 2 + h 2 r h 2 o 2 + h 

k -ll 

k_ 14 = 1.0 X 10 14 exp(-9800/RT) 

Recommendation by 
Richard S. Brokaw, 
NASA Lewis Research 
Center (1967) 

b 12 

k 12 

h 2 o 2 + M r OH + OH + M 

k -12 

k 12^N 2 = 1,7 X 1q17 exp(-46 300/RT) 

Baldwin and Brattan 
(ref. 17) 

13 

k 13 

h 2 + o 2 r OH + H 

k -13 

k 13 = 2.5 X 10 12 exp(-39 000/RT) 

Ripley and Gardiner (ref. 18) 


a All reactions have units of cm3/mole-sec except reactions 5 to 8, 10, and 12 have units of cm 6 / m ole 2 -sec; 
activation energies are expressed in units of cal/mole (1 cal/mole = 4.184 joules/mole). 
b See table 2 for relative third-body efficiencies. 
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TABLE 2.- RELATIVE THIRD-BODY EFFICIENCIES 



Source 

>Kaskan and Browne (ref. 15) 

J 

Wray (ref. 20) 

Bortner (ref. 21) 

to have an efficiency of unity) 

Kaskan and Browne (ref. 15) 

Baldwin and Brooks (ref. 22) 

*\ 

> Lewis and von Elbe (ref. 23) 

J 

Hoare, Pr other oe, and Walsh (ref. 24) 

\ Baldwin and Brattan (ref. 17) 

J 


III 
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TABLE 3.- TEST-GAS MIXTURE COMPOSITIONS AT 
ENTRANCE OF COMBUSTION CHAMBER 


Species 

Mass fractions for 


Mixture R 

Mixture A 
(a) 

Mixture B 
(b) 

N 2 

0.7568 

0.6191 

0.4794 

o 2 

.2162 

.2177 

.2192 

H 2 

.0270 

.0272 

.0274 

co 2 

0 

.0748 

.1507 

H 2 o 

0 

.0612 

.1233 


a Mixture A is based on a tunnel stagnation temperature of 


1390 K. 

^Mixture B is based on a tunnel stagnation temperature of 


2200 K, 



TABLE 4.- INITIAL CONDITIONS FOR CALCULATIONS 


Test 

case 

Mixture 

(a) 

Initial 

temperature, 

°K 

Pressure, 

atm 

1 

R 

1250 

1 

2 

A 





3 

B 




4 

R 



0.5 

5 

A 





6 

B 



T 

7 

R 



2 

8 

A 





9 

B 



f 

10 

R 

1500 

1 

11 

A 





12 

B 





13 

R 

1150 



14 

A 





15 

B 






a See table 3 for composition of mixtures R, A, and B. 



TABLE 5.- INITIAL MOLE FRACTIONS a OF O, H, AND OH 


Free radicals 


Initial mole fractions 2 

for - 


T 0 = 1150 K; 
p = 1 atm 

T 0 = 1250 K; 
p = 0.5 atm 

T 0 = 1250 K; 
p = 1 atm 

T 0 = 1250 K; 
p = 2 atm 

T 0 = 1500 K; 
p = 1 atm 

Mixture R: 






O 

2.67 X 10 " 9 

3.13 x 10-8 

2.21 X 10-8 

1.56 X 10-8 

1.28 X 10 " 6 

H 

2.12 X 10 ~ 8 

1.94 X 10 - 7 

1.38 X lO " 7 

9.83 X lO " 8 

5.04 X lO " 6 

OH ... . 

0 

0 

0 

0 

0 

Mixture A: 






0 

2.67 X 10 " 9 

3.13 x 10-8 

2.21 x 10-8 

1.56 x lO " 8 

1.28 x lO " 6 

H 

2.12 X 10-8 

1.94 x 10 - 7 

1.38 x 10 - 7 

9.83 x lO " 8 

5.04 X lO " 6 

OH ... . 

1.08 X 10 ~ 6 

4.96 X 10-6 

4.16 X lO " 6 

3.54 X 10-6 

5.59 X lO " 5 

Mixture B: 






0 

2.67 X 10 " 9 

3.13 X 10-8 

2.21 X 10-8 

1.56 X lO " 8 

1.28 X lO ” 6 

H 

2.12 X 10“8 

1.94 X 10 " 7 

1.38 X lO " 7 

9.83 X lO " 8 

5.04 X lO " 6 

OH ... . 

1.52 X 10-6 

7.04 X IQ -6 

5.89 X lO " 6 

4.97 X lO ' 6 

7.90 x 10 " 5 


a Based on equilibrium amounts in the two streams (hydrogen- and oxygen-containing 
streams) at the initial temperature T 0 and pressure p in the combustor prior to 
mixing. 



TABLE 6.- INDUCTION TIMES FOR VARIOUS INITIAL TEMPERATURES, 
PRESSURES, AND AMOUNTS OF C0 2 AND H 2 0 




Induction time, jusec, 

for — 


Mixture 

T 0 = 1150 K; 

T 0 = 1250 K; 

T 0 = 1250 K; 

T 0 = 1250 K; 

T 0 = 1500 K; 


p = 1 atm 

p = 0.5 atm 

p = 1 atm 

p = 2 atm 

p - 1 atm 



No initial O, H, and OH 



R 

62.6 

67.3 

34.5 

18.8 

13.0 

A 

76.4 

69.0 

36.8 

22.6 

13.5 

B 

102.2 

71.0 

40.0 

29.7 

13.2 



Initial amounts a of O, H, and OH 


R 

62.3 

65.9 

33.9 

18.5 

12.2 

A 

68.7 

57.0 

30.8 

19.6 

9.7 

B 

90.8 

56.8 

32.7 

25.1 

9.2 


a Based on equilibrium amounts in the two streams (hydrogen- and oxygen- 
containing streams) at the initial temperature and pressure in the combustor prior 
to mixing. 
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TABLE 7.- REACTION TIMES AND EQUILIBRIUM TEMPERATURES FOR VARIOUS 
INITIAL TEMPERATURES, PRESSURES, AND AMOUNTS OF CO 2 AND H 2 0 

(a) Reaction times 





r r , jusec, for - 



Mixture 

T 0 = 1150 K; 

T 0 = 1250 K; 

T 0 = 1250 K; 

T 0 = 1250 K; 

T 0 = 1500 K; 


p = 1 atm 

p = 0.5 atm 

p = 1 atm 

p = 2 atm 

p = 1 atm 

R 

65.4 

177.7 

56.5 

18.3 

43.0 

A 

77.6 

200 

67.2 

22.4 

47.5 

B 

91.8 

226 

73.0 

26.5 

54.8 


(b) Equilibrium temperatures 





T e , K, for - 



Mixture 

T 0 = 1150 K; 

T 0 = 1250 K; 

T 0 = 1250 K; 

T 0 = 1250 K; 

T 0 = 1500 K; 


p = 1 atm 

p = 0.5 atm 

p = 1 atm 

p = 2 atm 

p = 1 atm 

R 

2696 

2682 

2730 

2778 

2816 

A 

2588 

2582 

2625 

2666 

2712 

B 

2513 

2513 

2550 

2588 

2644 
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Species mass fractions 


Equilibrium 
values : 



2 468 2 468 2 468 

1 10 10 10 


Time, t, ysec 

Figure 2.- Time variation of species mass fractions and equilibrium levels for case 2 (T 0 = 1250 K, p = 1 atm, mixture A). 
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Species mass fractions 





Species mass fractions , 














Relative reduction in induction time, -At/t 


1 


10 


-1 


10 


-2 
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p(atm) Mixtures: 
R A B 
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0.5 □ a ■ 

2.0 O ❖ ♦ 


2 h 


10 


>3 


10 


-8 




4 


8.24 
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Sum of initial mole fractions of 0, H, and OH 


4 


I 


8 


10 
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Figure 18.- Relative reduction in induction time as a function of the sum of initial mole fractions of 0, H f and OH. (Flagged symbols 
are for a set of initial mole fractions taken to be less than equilibrium amount.) 
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